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(§j) Frequenzanalyseverfahren 

(g) Ein Frequenzanalyseverfahren umfaSt die Bewertung ei- 
nes zeitlichen, in Form von diskreten Abtastwerten vorlie- 
genden Eingangssignals mit einer Fensterfunktion. Anschlie- 
Send wird das gefensterte Eingangssignal zur Erzeugung 
eines Satzes von Koeffizienten Fourier-transformiert. 
Um ein solches Verfahren so weiterzubilden, daS die 
Eigenschaften des menschiichen Oh res nicht nur hinsicht- 
lich der spektralen Abbildung in den Frequenzbereich, 
sondern auch hinsichtlich der Auflosung im Zeitbereich 
nachgebildet werden, wird ein Block des Eingangssignals 
mit einem Satz von unterschiedlichen Fensterfunktionen zur 
Erzeugung eines Satzes von mit den jeweiiigen Fensterfunk- 
tionen gewichteten Blocken von Abtastwerten bewertet, 
deren Fourier-Transform ierte unterschiedliche Bandbreiten 
haben, bevor jeder der zeitgleich erzeugten Blocke von 
Abtastwerten einer eigenen Fourier-Transformation derart 
* unterzogen wird, daS fur jede Fensterfunktion zumindest 
J jeweils ein Koeffizient berechnet wird, dem die Bandbreite 
der Fourier-Transformierten dieser Fensterfunktion zugeord- 
net ist und daft die Koeffizienten so gewahlt werden, daS die 
ihnen zugeordneten Frequenzbander im wesentlichen anein- 
ander angrenzen. 
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Beschreibung 

Die vorliegende Erfindung betrifft ein Frequenzanalyseverfahren, bei dem aus einem zeitdiskret, in Form von 
BIdcken von Abtastwerten vorliegenden Eingangssignal mittels einer Fensterfunktion uberlappend aufeinander- 
5 folgende BIdcke von Abtastwerten entnommen werden, woraufhin die so mit den Fensterfunktionen bewerteten 
BIdcke zur Erzeugung eines Satzes von Koeffizienten fourier-transformiert werden, nach dem Oberbegriff des 
Patentanspruchs 1. 

Es sind hauptsachlich zwei verschiedene Frequenzanalyseverfahren bekannt, mittels derer ein zeitiiches, in 
Form von Abtastwerten vorliegendes Eingangssignal in eine Mehrzahl von Koeffizienten, d. h. in den Spektral- 

io bereich transformiert werden kann. Diese Verfahren lassen sich in sogenannte Teilbandverfahren (vgl. bspw. 
DE 40 1 1 577 Al) und Transformationsverfahren (vgl. bspw. US 51 42 656) kiassifizieren, wenngleich Transfor- 
mations- und Filterbankanalysen mathematisch zun&chst eng miteinander verwandt sind. 

FUr die Audiosignalanalyse und/oder -codierung ist es zweckma&ig und vorteilhaft, die Frequenzanalyse den 
Eigenschaften des menschlichen Ohres als Schallempfanger anzupassen. Eine wesentliche Erkenntnis der Psy- 

15 choakustik ist hierbei, daB das Ohr Schallereignisse, die in eine bestimmte Bandbreite, die sogenannte Frequenz- 
gruppenbreite fallen, gemeinsam bewertet Schallereignisse auBerhalb einer solchen Frequenzgruppenbreite 
werden getrennt wahrgenommen. 

Die Frequenzgruppenbreite ist — im Gegensatz zu rein technischen Frequenzanalysen — nicht konstant; mit 
zunehmender Mittenfrequenz wird auch ihre Bandbreite groBer. Der Zusammenhang zwischen Frequenzgrup- 

20 penbreite und Mittenfrequenz ist daruberhinaus, wie in Fig. 2 gezeigt, nichtlinear. 

Dabei verwenden die sogenannten Teilbandverfahren zur Spektralanalyse Filterbanke, die eine spektrale 
Zerlegung entsprechend den Frequenzgruppen vornehmen. Oblicherweise verwendet man hierzu Kaskadenan- 
ordnungen aus Halbbandfiltern. Durch die fortgesetzte Halbierung des Basisbandes mittels einer derartigen 
Halbbandfilterkaskade kann nur eine grobe Annaherung an die Aufteilung in Frequenzgruppen erreicht werden. 

25 Die Zeitaufldsung von Bandfiltern ist umgekehrt proportional zu ihrer jeweiligen Bandbreite; sie wird also mit 
zunehmender Bandbreite feiner und ist damit gut den Eigenschaften des Gehors angepaBt. Halbband-Filterkas- 
kaden erzeugen aber, bedingt durch das Prinzip der fortgesetzen Halbierung des Basisbandes, relativ groBe 
Abweichungen der Frequenzaufteilung von der Vorgabe durch die Frequenzgruppenaufteilung, die Zwicker fur 
das Gehdr gefunden hat. 

30 FFT-Analysen kdnnen zwar bei entsprechend feinem Frequenzraster eine genaue Aufteilung in Frequenz- 
gruppen gewahrleisten, die Zeitaufldsung ist dann jedoch zu grob und insbesondere konstant fur alle Frequenz- 
gruppen, was wiederum auch nicht den Eigenschaften des menschlichen Schallrezeptors entsprechen kann, wie 
die Wahrnehmbarkeit von Vorechostdrungen beweist 

Prinzipiell ist die Zeitaufldsung durch die Rekonstruktion einer Impulsfolge mit Hilfe eines idealen Tiefpasses 

35 gegeben. Die aus der Impulsantwort eines Filters ablesbare Zeitaufldsung wird als "resolution span" bezeichnet 
Sie ist in Fig, 1 wiedergegeben. 

Wie in Fig. 1 gezeigt ist, ist die Erzeugung der auf f g bandbegrenzten Funktion f(t) durch Superposition von 
mit dem jeweiligen Abtastwert bewerteten sinx/x-Funktionen a n (t) mdglich. Die sinx/x-Funktionen (Spaltfunk- 
tionen) stellen die Impulsantwort eines idealen Tiefpasses mit der Grenzfrequenz f g dar. Die Zeitaufldsung des 

40 idealen Tiefpasses ist T. 

Das Prinzip wird fiir ideale Interpolatoren verwendet Dazu werden sinx/x-Funktionen mit dem jeweiligen 
Abtastwert bewertet und urn einen Zeitabschnitt T verschoben, der sich aus dem Abstand von Maximum zum 
ersten Nulldurchgang der sinx/x- Funktion ergibt Die Superposition aller dieser sinx/x-Funktionen ergibt die 
interpolierte Zeitfunktion. Der Verlauf der auf diese Weise rekonstruierten kontinuierlichen Zeitfunktion stimmt 

45 an den Nulldurchgangen der Spaltfunktionen mit den Abtastwerten exakt uberein. Diese Bedingung ist als 
Interpolationsbedingung bekannt 

Als Zeitaufldsung wird eben dieser Abstand von Maximum zum ersten Nulldurchgang bzw. der Abstand der 
MuIldurchgSnge der sinx/x-Funktion bezeichnet Es ist leicht einzusehen, daB die Zeitaufldsung gleich der 
Abtastperiode ist, wenn als Interpolationsfilter ein TiefpaB mit einer Grenzfrequenz f g gleich der halben 

so Abtastfrequenz f a , also ein AllpaB gewahlt wird. Das Produkt aus Bandbreite und Zeitaufldsung betr&gt in 
diesem Fall B«T — f»- Vih — 0,5. Bei einer bestimmten Bandbreite kann die damit gegebene Zeitaufldsung nicht 
unterschritten werden. Dies wird in Analogic zu einem ahnlichen, aus der Physik bekannten Sachverhalt als 
Heisenberg-Grenze bezeichnet Teilband- und/oder Transformationsverfahren zur spektralen Zerlegung finden 
sich in fast alien heute bekannten Codern. Betde sind fur eine gehdrgerechte Representation im Frequenz- und 

55 Zeitbereich ungeeignet Speziell muB ein MeBverfahren, das beansprucht, die QualitSt heutiger Codierverfahren 
beurteilen zu kdnnen, sicher bereits hinsichtlich der Frequenzanalyse genauer arbeiten; erst dann ist die Genau- 
igkeit, mit der hieruber die Mithdrschwelien bestimmt werden hoher, als in den zu beurteilenden Coder- Deco- 
der- Kombinationen. 

Die Nachteile herkdmmlicher, heute fast ausschlieBlich angewendeter Verfahren zur gehdrrichtigen Fre- 
60 quenzanaiyse seien in zusammengefaBter Form verdeutlicht: 

— Teilbandverfahren arbeiten zwar an der Heisenberg-Grenze; die Frequenzgruppenbreiten werden 
jedoch nur grob getroffea 

— Transformationsverfahren sind zwar in der Lage, die Frequenzgruppenbreiten im Prinzip beliebig genau 
65 nachzubilden, die Zeitaufldsung ist jedoch bei weitem zu grob, daruber hinaus konstant und somit nicht 

gehdradaquat 

— Hybride Frequenzanalysen nahern sich der geforderten Zeit-/Frequenzaufldsung des Gehors zwar 
besser, eine exakte Nachbildung ist aber auch hiermit prinzipiell nicht moglich. 
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Der Erfindung liegt die Aufgabe zugrunde, ein Frequenzanlayseverfahren zu schaffen, das die Eigenschaften 
des menschlichen Ohres nicht nur im Frequenzsondern gleichzeitig auch im Zeitbereich nachbildet Diese 
Aufgabe wird durch ein Frequenzanaiyseverfahren gemafl Anspruch 1 gelost 

GemaB der Erfindung werden bei einem Frequenzanaiyseverfahren der eingangs genannten Art aus einem 
zeitlich in Form von Abtastwerten vorliegenden Eingangssignai mittels Fensterfunktionen Qberlappend aufein- 5 
anderfolgende Blocke von Abtastwerten entnommen, worauf die so mit den Fensterfunktionen bewerteten 
Blocke zur Erzeugung eines Satzes von Koeffizienten Fourier-transformiert werden, wobei aber jeweils ein 
soicher Block des Eingangssignales mit einem Satz von unterschiediichen, symmetrischen Fensterfunktionen 
gleicher Lange zur Erzeugung eines Satzes von mit den jeweiligen Fensterfunktionen gewichteten Bldcken von 
Abtastwerten bewertet wird Die Fourier-Transformierten der Fensterfunktionen besitzen, bedingt durch die 10 
unterschiedlichen Fensterfunktionen, auch unterschiedliche Bandbreiten. Zu jedem der zeitgleich erzeugten 
Blocke von gefensterten Abtastwerten laBtsich iiber eine eigene Fourier-Transformation mindestens ein Koeffi- 
zient berechnen, dem die Bandbreite der Fourier-Transformierten dieser Fensterfunktion zugeordnet ist, wobei 
die Koeffizienten so gewahlt werden, daB die ihnen zugeordneten Frequenzbander, bezogen auf ihre Bandbreite, 
sich nur geringffigig iiberiappen oder beabstandet sind Mit anderen Worten schlieBen die den Koeffizienten 15 
zugeordneten Frequenzbander im wesentlichen luckenlos und uberdeckungsfrei aneinander an. 

Ein bevorzugtes Ausfuhrungsbeispiel des erfindungsgemaBen Frequenzanalyseverfahrens wird nachfolgend 
unter Bezugnahme auf die beiltegenden Zeichnungen naher erlautert Es zeigen: 

Fig. 1 eine Darstellung zur Verdeutlichung, daB die Erzeugung einer auf Fg bandbegrenzten Funktion f(t) 
durch Oberlagerung von mit dem jeweiligen Abtastwert bewerteten sinx/x-Funktionen moglich ist; 20 

Fig, 2 den Zusammenhang zwischen der Frequenzgruppe AFg in Abhangigkeit von der Mittenfrequenz f m ; 

Fig. 3 den Zusammenhang zwischen einer Frequenzachse in Hertz als Ordinate und einer Frequenzachse in 
BARK ais Abszisse; 

Fig. 4 und 5 Frequenz- und Zeitdarstellungen zur Verdeutlichung des Ahnlichkeitssatzes der Fourier-Trans- 
formation zur Verdeutlichung, daB zu einer Verbreiterung der Bandbreite eines Spektralkoeffizienten eine 25 
Verschmalerung der Fensterfunktion im Zeitbereich ndtig ist; 

Fig. 6 ein Beispiel fur eine Fourier-Transformierte einer Fensterfunktion, welche die Energieerhaltung ge- 
w&hrleistet; 

Fig. 7 ein Beispiel fur eine Fourier-Transformierte, die eine feinere Zeitauflosung, jedoch nur eine reduzierte 
Frequenzaufldsung besitzt; 30 

Fig. 8 ein FluBdiagramm eines BARK-Generators zur Ermittlung von Fensterfunktionen im Zeitbereich, die 
frequenzgruppenbreite Spektralkoeffizienten erzeugen; 

Fig. 9 eine Darstellung zur Beurteilung eines BARK-Generatorentwurf es; 

Fig. 1 0 eine Darstellung des Iterationsergebnisses des BARK-Generators in Abhangigkeit von der Rasterwei- 
te; 35 

Fig. 1 1 eine Darstellung des Quotienten aus der Anzahl unterschiedlicher Frequenzgruppen und der Gesamt- 
zahl der Bander in Abhangigkeit von der Rasterzahl; 

Fig. 12 ein Blockdiagramm einer Vorrichtung zur Durchfiihrung der BARK-Transformation auf der Grundla- 
ge einer 51 2-Punkte-FFT; 

Fig. 1 3 ein rekursives Filter erster Ordnung zur Bildung des Goertzel-Algorithmus; 40 

Fig. 14 ein rekursives Filter zweiter Ordnung zur Bildung eines Goertzel-Algorithmus zweiten Grades; 

Fig. 1 5 eine transponierte Struktur eines Filters zur Bildung eines Goertzel-Algorithmus zweiten Grades; 

Fig. 16 eine bevorzugte Ausfuhrungsform eines rekursiven Filters zur Bildung des Goertzel-Algorithmus 
zweiten Grades als Vorrichtung zur Berechnung eines Koeffizienten der BARK-Transformation; 

Fig. 17 eine Darstellung zum Vergleich des Rechenaufwandes fur einen komplexen Spektralkoeffizienten bei 45 
Verwendung einer schnellen Fourier-Transformation gegenuber dem Goertzel-Algorithmus, aufgetragen iiber 
verschiedene Transformationslangen N; 

Fig. 18 eine Darstellung der Anzahl GA der Goertzel-Algorithmen, die innerhalb des Rechenzeitbedarfs fur 
eine schnell Fourier-Transformation auf einem Signalprozessor gerechnet werden konnen, aufgetragen iiber 
verschiedene Transformationslangen N; 50 

Fig. 19 eine Darstellung der Betragsubertragungsfunktion der Einzelfilter einer BARK-Transformation; und 

Fig. 20 eine Darstellung des Frequenzganges als Summe der Amplitudenubertragungsfunktionen. 

Transformationsverfahren Iassen sich offensichtlich mit deutlich weniger Aufwand als Filterbanke realisieren. 
Auf der anderen Seite sind sie aber nur dann mit Erfolg einsetzbar, wenn auf der spektralen Achse (der 
Frequenzachse in Hertz) lediglich aquidistant angeordnete Reprasentanten bendtigt werden. Eine nicht aquidist- 55 
ante Anordnung auf der spektralen Achse iiber eine Gruppierung herzustelien, fiihrt zu einer Realisierung, die 
von vorne herein nicht gehorrichtig ist, weil sie nicht entlang der Heisenberg-Grenze arbeitet 

Der Zusammenhang zwischen einer Frequenzachse in Hertz und einer Frequenzachse in BARK ist aber — 
jedenfalls fiir niedrige Frequenzen — naherungsweise linear. Fig. 2 [Zwicker, E.; Feldtkeiler, R.: Das Ohr als 
Nachrichtenempf anger, 2. Aufl, Hirzel Veriag, Stuttgart 1967] laBt dies bereits erkennen. Deutlicher wird der eo 
lineare Zusammenhang, wenn man die Frequenz in Hertz uber der Tonheit in BARK auftragt, wie Zwicker dies 
in Fig. 3 getan hat [Zwicker, E.; Feldtkeiler, R-: Das Ohr als Nachrichtenempfanger, 2. AufL, Hirzel Veriag, 
Suttgart 1967]. Dann ist erkennbar, daB der Zusammenhang zwischen einer Frequenzachse f in Hertz (Ordinate) 
und einer Frequenzachse z in BARK (Abszisse) zumindest im unteren Frequenzbereich bis etwa 5 BARK linear 
ist; das bedeutet, das hier die Teilbander gleich breit sind. Im Bereich von 0 Hz bis etwa 500 Hz kann also in guter 65 
Naherung von einem linearen Zusammenhang gesprochen werdea In diesem Bereich laBt sich daher fur jede 
Mittenfrequenz in etwa die gleiche Frequenzgruppenbreite von 100 Hz angeben. Hier wird zunachst eine 
Fourier-Transformation vom Zeitbereich in den Frequenzbereich eingesetzt Sie kann hier ohne die iibliche 
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Gruppierung von Spektralwerten, d h. mit einem idealen Produkt von Bandbreite und Zeitaufldsung verwendet 
werden. Man gewinnt dabei den Vorteil, zumindest die ersten funf Spektralkoeffizienten gemeinsam in einem 
Ablauf (z. B. in einer schnellen Diskreten Fourier-Transformation) berechnen zu konnen. Oblicherweise ergeben 
sich bei der Abbildung einer Frequenzachse von 0 Hz bis etwa 22 kHz auf eine BARK-Achse etwa 25 Spektralli- 
5 nien. Sie kdnnen als Reprasentanten fur die frequenzgruppenbreiten Bander angesehen werden. Die Berech- 
nung der ersten funf Spektrallinien mittels einer Fourier-Transformation hat mit einer Schrittweite auf der 
Frequenzachse zu erfolgen, die der in diesem Bereich annahernd konstanten Frequenzgruppenbreite, also etwa 
100 Hz, entspricht Bei einer Abtastfrequenz von beispielsweise 48 kHz muB demnach mit einer 480-Punkte 
Fourier-Transformation gerechnet werden. Von den 240 Spektrallinien kdnnen dann schon die ersten 5 verwen- 

io det werden, die restlichen 235 Spektrallinien mOssen verworfen werden, da sie zu hdheren Frequenzen hin keine 
Frequenzgruppenbreite mehr reprasentieren (ihre Bandbreiten sind dafur zu schmal). Fur den geforderten 
Entwurf entlang der Heisenberg-Grenze sind wegen der Reziprozitat von Bandbreite und Zeitaufldsung B -T « 
0,5 zur Erzielung breiterer Bander entsprechend schmalere Fenster im Zeitbereich fur weitere, getrennt zu 
berechnende Fourier-Transformationen zu wahlen; diese ergeben dann — wie gewQnscht — auch eine bessere 

15 Zeitaufldsung. 

Der Ahnlichkeitssatz der Fourier-Transformation besagt, daB zu einer "Dehnung" der Bandbreite eines 
Spektralkoeffizienten eine "Stauchung" der Fensterfunktion im Zeitbereich ndtig ist Zeitaufldsung T und Band- 
breite B verhalten sich in reziproker Weise zueinander. Diese Reziprozitat wird mittels der in Fig. 4 und Fig. 5 
veranschaulichten Korrespondenzen der Fourier-Transformation gezeigt Zur Erzielung einer hdheren Band- 

20 breite muB demnach die Fensterfunktion im Zeitbereich gestaucht werden; seine Zeitaufldsung wird damit 
besser. Der dem erfindungsgemaBen Frequenzanalyseverfahren zugrunde liegende Gedanke ist es also, die 
Ordnung der Fourier-Transformation (die Anzahi der Koeffizienten) beizubehaltea gieichzeitig aber dafur zu 
sorgen, daB die hdheren Spektrallinien Reprasentanten fur breitere Bander darstellen. Auf diese Weise erreicht 
man auch gleiche Steilheiten fQr die Filterflanken. Dies ist die Voraussetzung fur ein konstantes GesamtQbertra- 

25 gungsverhalten fur den Betragsfrequenzgang, das fur den Einsatz in MeBverfahren, insbesondere aber auch in 
Codierverfahren, gefordert werden muB. 

Mit einer BARK-Transformation wird eine gehdrgerechte Transformation vom Zeitbereich in den Frequenz- 
bereich bezeichnet, wobei eine aquidistante Darstellung der Spektralkoeffizienten auf einer BARK-Achse (und 
nicht auf einer Hertz-Achse) erreicht wird. Die hier entwickelte BARK-Transformation unterscheidet sich von 

30 anderen (meist nur gehdrangepaBten) Abbildungen auf die BARK-Achse dadurch, daB sie entlang der Heisen- 
berg-Grenze und damit gehdrrichtig arbeitet 

Eine Fourier-Transformation, die grdBere -6 dB- Bandbreiten liefert, als es dem Abstand Af der Spektralli- 
nien entspricht, ist im herkdmmlichen Sinne naturlich unbrauchbar (Fig. 6). 
Normalerweise ist die Bandbreite Bn aller Spektralkoeffizienten der Fourier-Transformation durch die Ord- 

35 nung N und die Abtastfrequenz f a vorgeschrieben: B n =» f a /N. Die Betragsubertragungsfunktionen schneiden 
sich dann in |H(f)| - 03, d. h. an den -6 dB-Punkten (Fig. 7). Dennoch konnen von einer nach Fig. 6 modifizier- 
tcn Fourier-Transformation nur die Linien fur die BARK-Transformation verwendet werden, deren Bandbreiten 
bzw. Mittenfrequenzen in das von Zwicker fur die Frequenzgruppenbreiten gegebene Schema passea Durch die 
— zumindest innerhaib eines Rasters — freie Wahl der Bandbreite jedes Koeffizienten eriibrigt sich eine 

40 Gruppierung von einigen Koeffizienten zu einer Frequenzgruppe mit all ihren Nachteilen, wie sie bei herkdmm- 
lichen gehdrangepaBten Transformationsverfahren ndtig ist Weil das Frequenzraster fur jede neu zu suchende 
Linie einer BARK-Transformation beibehaiten wird, lassen sich nun Bander aneinanderreihen, deren Breite in 
beliebig guter Naherung der angestrebten Abbildung auf die BARK-Achse entspricht Ein Entwurf entlang der 
Heisenberg-Grenze ist dabei gewahrleistet Die Staffelung ist deshalb "luckenios" mdglich, weil Bandbreiten 

45 gewahlt werden, die genau in das durch die Lange der Fourier-Transformation gegebene Af-Raster passen. Die 
Summe aller AmplitudenObertragungsfunktionen bleibt dabei dennoch konstant: die Flankensteilheit der bei der 
Staffelung entstehenden Bandfilter ist nur von der Anzahi der Filterkoeffizienten (der Fensterlange), nicht aber 
von der Bandbreite abhangig. Diese Anzahi ist aber konstant, sie entspricht normalerweise (d. h. ohne Polypha- 
senaddition) der L3nge der Fourier-Transformation. Weil alle Flanken die gleiche Steilheit aufweisen, treten 

50 unerwflnschte Schwankungen im Gesamtubertragungsverhalten (Summenfrequenzgang) hier nicht auf. 

Die in Form von Abtastwerten vorliegenden Fensterfunktionen mussen nicht notwendigerweise die gleiche 
Anzahi von Abtastwerten aufweisen, die der Lange der auf die Fensterung folgenden Fourier-Transformation 
entspricht Vielmehr kann durch eine Verlangerung der Fensterfunktion uber die Transformationslange hinaus 
eine Verbesserung der Flankensteilheit der Einzelfilter erreicht werden. Das Additionstheorem der Fourier- 

55 Transformation erlaubt es, auch iangere Sequenzen als durch die Transformationslange gegeben in den Fre- 
quenzbereich zu transformieren. Hierzu wird eine Anzahi von aufeinanderfolgenden Teilbldcken aus dem 
bereits gefensterten langeren Eingangsblock von Abtastwerten gebildet, nach dem Prinzip der Superposition 
addiert und gemeinsam in den Frequenzbereich transformiert Dieses Prinzip ist auch als Polyphasen-Filterung 
bekannt[Chrochiere, R.E.; Rabiner L: Multirate Digital Signal Processing, Prentice- Hall, Englewood Cliffs, New 

60 Jersey, 1983], Durch die Polyphasenfilterung wird insbesondere die Zeitaufldsung nicht berilhrt, da die dem 
Filterentwurf zugrunde liegende Spaltfunktion die gleiche bleibt 

Die Erfindung liefert ferner ein als BARK-Generator bezeichnetes Iterationsverfahren, das nunmehr erlautert 
wird. Die Staffelung der Bander wird nach Vorgabe einiger Parameter einem Rechenautomaten Qberlassen. 
Diese Staffelung liefert fur jede Frequenzgruppe einen Satz von Filterkoeffizienten fur das Zeitfenster einer 

65 schnellen Fourier-Transformation. Es hat sich als gunstig erwiesen, die Staffelung ab der halben Abtastfrequenz 
(d. h. der hdchsten vorkommenden Signalfrequenz) zu tieferen Frequenzen hin und nicht umgekehrt vorzuneh- 
men. (Eine Staffelung ab der niedrigsten im Signal vorkommenden Frequenz ist naturlich auch mdglich). Die 
Staffelung gewahrleistet nicht von vorne herein, daB am Ende die untere Grenzfrequenz der untersten Fre- 
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quenzgruppe mit 0 Hz zusammenfallt Brauchbar sind aber auch Iterationen, die untere Grenzfrequenzen fur das 
unterste Band bis hinauf zu etwa 50 Hz liefern. Ein moglicherweise im Signal enthaltener, normaierweise 
unhorbarer Gleichanteil wird dadurch in die Messung stark gedampft eingehen. Audiosignale reichen nicht bis 
zu einer Frequenz bis 0 Hertz hinab. Im Verlauf der Iteration werden also weitere Bander unter. die bereits 
ermittelten gestaffeit Die Bandgrenzen der frequenzgruppenbreiten Bander, aber auch deren Mittenfrequenzen 5 
fallen dabei auf ein vorzugebendes Raster. 

Bei der Suche nach den Koef fizientensatzen fur weitere Fensterfunktionen wird nach dem in Fig. 8 dargestell- 
ten FluBdiagramm vorgegangen (CBW = Frequenzgruppenbreite nach Zwicker). 

— Die obere Bandgrenze ist bekannt Ste ergibt sich aus der unteren Bandgrenze des bereits als passend 10 
ermittelten Zeitfensters fur die vorherige, hohere Frequenzgruppe, wobei 6 dB-Bandbreiten zugrunde 
gelegt werden; bei Iterationsbeginn wird die obere Grenzfrequenz auf die halbe Abtastfrequenz als Start- 
wert gesetzt 

— Die Mittenfrequenz des als nachstes zu suchenden Bandes wird nun solange auf dem Raster verschoben, 
das sich durch die Abtastfrequenz und die Ordnung der Fourier-Transformation ergibt,. bis aus dem 15 
Ensemble von moglichen Bandbreiten eine gefunden wird, die der Frequenzgruppenbreite von Zwicker am 
nachsten kommt Ist das verwendete Frequenzraster nur fein genug, so wird auch die Genauigkeit, mit der 
eine Frequenzgruppenbreite getroffen wird, beliebig gut sein. 

— Aus der gefundenen Bandbreite ergibt sich der Koeffizientensatz fur eine der Fensterfunktionen und 
damit die Zeitauflosung dieses Fensters gemaB der Heisenberg-Grenze, d h. die Anforderung im Frequenz- 20 
bereich bestimmt dabei den Verlauf der Fensterfunktion im Zeitbereich. Hier wird das Filterentwurfsver- 
fahren nach Tietze, U.; Schenk, Ch.: Halbleiterschaltungstechnik, 9. Auflage, Springer- Verlag, Berlin, 1989, 
eingesetzt 

— Die untere Grenzfrequenz des gefundenen Bandes .bildet die obere Grenzfrequenz des als nachstes zu 
ermittelnden Bandes. Wird im Verlauf der Iteration die Mittenfrequenz negativ, so ist das Abbruchkriterium 25 
erreicht 

Es ist uberraschend, daB bei der Untereinanderstaffelung der frequenzgruppenbreiten Bander weit weniger 
als 25 unterschiedliche Zeitfenster notig sind. Dies ist einerseits auf die nahezu konstante Frequenzgruppenbrei- 
te zu tiefen Frequenzen bis etwa 5 BARK zuruckzufiihren (vgl. Fig, 3); andererseits ist aber auch entscheidend, 30 
daB durch die Vorgabe eines Frequenzrasters eben nur diskrete Mittenfrequenzen bzw. Bandgrenzen erlaubt 
sind. Wie Fig. 9 zeigt, konnen bereits bei einer Diskretisierung auf 512 Frequenzen bis zur halben Abtastfre- 
quenz die 25 frequenzgruppenbreiten Bander ausreichend genau mit nur 14 verschiedenen Fensterfunktionen 
nachgebildet werden: Die Beurteilung des BARK-Generator- Entwurfs ist sofort graphisch am Bildschirm mog- 
lich. Die Abweichung der Bandbreiten von den Frequenzgruppenbreiten nach Zwicker ist bereits bei einem 35 
relativ groben Frequenzraster von 512 diskreten Frequenzen ausreichend klein (Beispiel in Fig. 9 fur 44,1 kHz 
Abtastfrequenz). 

Der BARK-Generator ieistet die folgenden Aufgaben: 

— Der Entwurf ist prinzipiell fur jede Abtastfrequenz moglich, speziell aber fur die gangigen Abtastfre- 40 
quenzen 44,1 kHz, 48 kHz und 32 kHz; eine Abtastfrequenz von 44,1 kHz ist fur Compact- Disc Musikmate- 

rial ublich; Abtastf requenzen von 32 kHz und 48 kHz findet man bei DAT-Medien (Digital Audio Tape). 

— Untereinanderstaffelung von frequenzgruppenbreiten Bandera. Die Genauigkeit, mit der der Entwurf 
den Frequenzgruppenbreiten von Zwicker entspricht, laBt sich durch die Vorgabe eines Rasters steuern. 

— Zu Versuchszwecken lassen sich auch Bandbreiten fur die Teiifilter einstellen, die Bruchteilen von einem 45 
BARK entsprechen; dies eroffnet die Moglichkeit, das Tonhdhenauflosungsvermogen des menschlichen 
Ohres von etwa 1/27 BARK nachzubilden. 

— Berechnung und Ausgabe der Satze von Filterkoeffizienten fur die bendtigten Fenster unter BerGcksich- 
tigung der Polyphasenaddition in verschiedenen Ausgabeformaten. Eine der Ausgabemoglichkeiten ist eine 
Konstantentabelle, die sofort als Quelltext in ein Programm zur BARK-Transformation eingebunden 50 
werden kann. Die Bezeichnung BARK-Generator ist damit gerechtfertigt 

— Bandbreite und Zeitauflosung der Fensterfunktionen ergeben (ahnlich wie bei der trivialen Bandfilter- 
Kaskade) immer einen Entwurf entlang der Heisenberg-Grenze. 

— Eine visuelle Beurteilung der Iterationsergebnisse ist sofort durch die eingebundene Graphikausgabe 
moglich. Eventuelle Modifikationen konnen interaktiv erfolgen. 55 

Wie Fig. 10 zeigt lauft die Iteration nicht fur beliebige Frequenzraster, sondern nur bei einigen ausgezeichne- 
ten Rasterweiten erfolgreich ab. Die Iteration kann immer dann als erfolgreich gelten, wenn die untere Grenz- 
frequenz des untersten Bandes im Bereich von etwa 0 bis 50 Hz Hegt Es ist giinstig, einen im Signal mdglicher- 
weise enthaltenen Gleichanteil, der ja normaierweise unhdrbar ist, auch in die Messung stark gedampft eingehen 60 
zu lassen. Der DurchlaBbereich der untersten Frequenzgruppe mufl daher nicht unbedingt bis 0 Hz hinabrei- 
chen. Die Gesamtzahl der Bander, in die der Freqiienzbereich von z. B. 0 bis 22,05 kHz unterteilt wird, n^hert sich 
mit zunehmender Rasterzahl (d. h, mit feiner werdendem Frequenzraster) der Frequenzgruppenzahl von 25. 
Diese Anzahl ist fur die in diesem Beispiel gewahlte Abtastfrequenz richtig. 

Fig. 10 zeigt auch, daB die Anzahl unterschiedlicher Fensterfunktionen dabei immer unterhalb der gesamten 65 
Anzahl von Frequenzgruppen bleibt Hierin Hegt ein wesentlicher Vorteil des BARK-Generators: Es sind — fur 
alle praktisch relevanten Falle — weniger Bewertungen des Eingangsvektors des Zeitsignals mit einer Fenster- 
funktion notig als durch die Transformation frequenzgruppenbreite Bander erzeugt werden. Bildet man den 

5 
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Quotienten aus der Anzahl unterschiedlicher Frequenzgruppen und der-Anzahi der Bander, so ergibt sich fur 
eine Rasterzahl von etwa 512 ein Minimum (Fig. 11). 

Hier arbeitet die BARK-Transformation am okonomischsten: Im Verhaltnis wird die geringste Anzahl unter- 
schiedlicher Bander bendtigt Noch kleinere Rasterzahlen sind zwar moglich, fuhren in diesem Beispiel aber zu 
5 Abweichungen von der Frequenzgruppenbreite nach Zwicker, die nicht mehr befriedigen. 

MaBnahmen zur Reduktion des Rechenaufwandes 

Der Rechenaufwand zur Ermittiung einer vollstandigen BARK-Transformation ist noch relativ groB, da fur 
io jedes unterschiedliche Fenster eine neue Diskrete- Fourier-Transformation gerechnet werden muB. Diese Vor- 
gehensweise geht mit einem meist nur begrenzt zur Verfflgung stehenden Potential an Rechenleistung offen- 
sichtlich sehr unokonomisch um, da die meisten der berechneten Fourier- Koeffizienten verworfen werden; nur 
sehr wenige Koeffizienten (in vielen Fallen sogar nur ein einziger) werden fiir die BARK-Transformation 
weiterverwendet 

15 Wie in Fig. 12 angegeben, kann eine BARK-Transformation, deren Gmndlage eine 512-Punkte FFT bildet, 
bereits mit nur 1 1 unterschiedlichen Teilfenstern berechnet werden. Bei der Berechnung der Koeffizienten fQr 
die 25 frequenzgruppenbreiten Teilfilter Ober eine schnelle Fourier-Transformation (FFT) mOssen aber 99,56% 
aller FFT- Koeffizienten verworfen werden, nur 25 aus 5632 berechneten Koeffizienten werden tatsachlich 
verwendet 

20 Fig. 12 verdeutlicht die Verschwendung von Rechenkapazitat in einem Beispiel mit 11 unterschiedlichen 
Fenstern und einer Rasterung auf 512 diskrete Frequenzen bis zur Abtastfrequenz: Fur die untersten 10 
Spektrallinien kann wegen des in diesem Frequenzbereich naherungsweise linearen Zusammenhangs zwischen 
Frequenzgruppenbreite und Mittenfrequenz noch ein und dasselbe Zeitfenster verwendet werden. Fiir die 
folgenden 3 Bander (11 bis 13) wird ein neues Zeitfenster und damit eine neue gemeinsame FFT bendtigt Dies 

25 gilt auch noch fiir die darauffolgenden 2 Frequenzgruppen; erst ab der 16ten Frequenzgruppe aufwarts muB 
jeweits ein neues Zeitfenster pro Frequenzgruppe eingesetzt werden. Bei der ersten FFT werden dadurch 502 
von 512 Koeffizienten verworfen, dann 509 von 512, 510 von 512 und schlieBIich 8 mal 511 von 512. Insgesamt 
werden von 1 1*512 - 5632 Koeffizienten nur 25 verwendet Es mussen also 5607/25 - 224 mal mehr komplexe 
Fourier-Koeffizienten berechnet werden, als dann far die BARK-Transformation gebraucht werden. Wie in 

30 Fig. 12 jedenfalls beispielhaft skizziert wird, ergabe die Berechnung der BARK-Transformation uber eine 
Anzahl von schnellen Fourier-Transformationen mit jeweils unterschiedlichen Zeitfenstern einen unverhaitnis- 
maBig hohen Rechenaufwand, wenn nicht eine effizientere Methode zur Berechnung einzelner Koeffizienten als 
die schnelle Fourier-Transformation eingesetzt wurde. 

Eine triviale Moglichkeit, nur einen Koeffizienten einer Fourier-Transformation zu berechnen, besteht darin, 

35 die Bestimmungsgleichung fQr eine diskrete Fourier-Transformation anzuwenden. Mehr Vorteile bietet aber der 
Einsatz des Goertzel-Algorithmus zu Berechnung einzelner Linien. Mehr Vorteile bietet aber der Einsatz des 
Goertzel-AIgorithmus zur Berechnung einzelner-Linien [Goertzel, G.: An Algorithm for the Evaluation of Finite 
Trigonometric Series, Amer. Math. Monthly 65 ( 1 968), S34— 35]. 

Ausgehend von der z-Transformation X(z) einer Sequenz x(n) von Abtastwerten der Lange N 

40 

n=-0 

45 

ergibt sich die Diskrete Fourier-Transformation (DFT) mit 

z - W k - ( e -j 2n/N ) k - e-i 2nk/N . 

50 Ober k werden diskrete Frequenzen ausgewahlt 

Mit der Substitution z n - W nk ~ e" j 2n " N kann die DFT auch als Polynom angegeben werden. Beispielsweise 
ergibt sich fOr N — 4 ein Polynom dritter Ordnung 

X(z) - xpj-z 3 + x(2)z 2 + x(l) z + x(0). 

55 

Ober das Horner-Schema Idst man die Multiplikationen mit Potenzen von z auf 

X(z) - «x(3)-z + x(2)).z + x(l))z + x(0). 

60 Unter Einfflhrung einer neuen Variablen y(n) kann daraus eine lineare Differenzengleichung angegeben 
werden 

y(n)-Z'y(n-l) + x(N-n). 

es Anfangsbedingung ist hierbei y(0) = 0. Das gesuchte Ergebnis fur die Fourier-Transformierte findet man fQr n 
- N 

X(z) - y(N). 

6 
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Die lineare Differenzengleichung kann auch als rekursives Filter erster Ordnung aufgefaBt werden. Die 
eingangs gemachte Substitution z n = W nk wird hierzu wieder ruckgangig gemacht 

y(n) s W k y(n-1) + x(N-n);y(0) = 0. 

5 

Die damit gegebene Anordnung ist in Fig. 13 wiedergegeben. 

Die Eingangssequenz muB allerdings in umgekehrter Reihenfolge vorliegen. 

AuBerdem wirkt sich nachteilig auf die Stabilitat und den Rechenaufwand filr das Filter aus, daB die Zustands- 
variable y(n-l) uber einen komplexen Koeffizienten W k auf den Filtereingang zuriickgefuhrt wird. 

Die Anzahl der Multiplikationen und Additionen ist zwar die gleiche wie bei der direkten Berechnung der io 
diskreten Fourier-Transformation, die Drehfaktoren miissen hier jedoch nicht explizit berechnet oder uber eine 
Drehfaktortabelle ("twiddle table") zur Verfugung gestellt werden; durch die rekursive Berechnungsweise 
werden sie vielmehr in jedem Durchlauf, ausgehend von einem Anfangswert, fflr einen neu hinzukommenden 
Abtastwert der Eingangssequenz x(n) weitergedreht 

Die Anzahl der bendtigten Multiplikationen laBt sich auf die Halfte reduzieren, wenn man die Euler'sche 15 
Gleichung einsetzt 

W k « e -J 2n/N = cos(27ik/N) - j sin(27ik/N) 

und auf eine indizierte Notation ubergeht 20 

y n = W k y n .i + XN-n « [cos(27tk/N) - jsin(27ik/N)] y n .| + xn-h. 

wobei y n -i = [cos(2nk/N) — j sin(2rck/N)]-y n .2 + xn-o+i- 

Nach einigen weiteren Umformungen ergibt sich schlieBlich fur Yn 25 

y n - [2 • cos(2ulc/N)] . y„_i + y n -2 = -xn-h + [cos(2u.k/N) + j sin (2uic/N)] ■ XN + i-n- 

Durch Vergleich mit der allgemeinen Differenzengleichung zweiter Ordnung 

y(k-t-2) + c, • y(k + 1) + Co • y(k) = b 2 • x(k+2) + bi . x(k+ 1) + b 0 * v(k) 



h u (t)=-h(t). 



30 



35 



kdnnen die Filterkoeffizienten fur ein rekursives Filter zweiten Grades in der ersten kanonischen Gruhdform 
angegeben werden, das sich wie in Fig. 1 4 gezeigt, darstellen laBt. 

ci » — 2-cos(2nk/N);co - 1 
b 2 - 0;bi - (27tk/N) + j sin (2nk/N); b 0 = 1. 

In dieser Form besitzt die Gleichung reelle Koeffizienten fur die Ruckfuhrung, allerdings einen komplexen 
Eingangskoeffizienten. 40 

Das rekursive Filter zweiten Grades in der zweiten kanonischen Struktur kann — wie in Fig. 15 gezeigt — 
noch transponiert werden. 

Da die Eingangsfolge rein reell ist, aber in umgekehrter Reihenfolge am Eingang anliegt, gilt fur den ungera- 
den Anteil der Eingangsfolge 



45 



Der ungerade Anteil h u (t) Iiefert aber nur den ungeraden Imaginarteil der Fourier-Transformation; speist man 
die Eingangsfolge in nicht umgekehrter Reihenfolge ein, muB der Imaginarteil der Fourier-Transformierten 
Iediglich negiert werden. 50 

Fig. 16 zeigt die endgiiltige Version des Goertzel-Algorithmus. Die Eingangsfolge x(n) liegt jetzt in normaler 
Reihenfolge am Eingang an. Der Ausgangswert muB erst bestimmt wenn der Eingangsvektor x(n) vollstandig 
abgearbeitet ist Zur Berechnung einer Spektrallinie ist es vorteilhaft, zunachst alle Werte der Eingangsfolge x(n) 
in das Filter "hineinlaufen" zu lassen und erst dann den Ausgangswert y(n) uber eine komplexe Multiplikation aus 
den Zustandsvariablen At und A2 zu berechnen. 55 

Mit diesem wenig bekannten Algorithmus ist es also moglich, einzelne Linien einer diskreten Fourier-Trans- 
formation zu berechnen; die dabei verwendete Struktur laBt sich auch als rekursives Filter auffassen. Ein 
besonderer Vorteil gegenuber der FFT ist, daB hier keine Drehfaktortabelle bendtigt wird AuBerdem treten 
innerhalb des Filters Iediglich rein reelle Additionen und Multiplikationen auf. Nur am Filterausgang wird eine 
komplexe Multiplikation bendtigt; zur Berechnung eines BARK- Koeffizienten muB diese auBerdem nur einmal eo 
berechnet werden. Der Aufwand zur Berechnung eines Koeffizienten uber den Goertzel-Algorithmus ist in der 
praktischen Realisierung auf einem modernen Einchip-Signalprozessor TMS320C30 geringer als der far eine 
komplette FFT. (Fig. 1 7 zeigt den Rechenaufwand fflr verschiedene Transformationsiangen N). 

Fig. 18 zeigt dazu, daB innerhalb der Rechenzeit die eine schnelle Fourier-Transformation bendtigt, in vielen 
Fallen sogar einige Goertzel-Algorithmen gerechnet werden konnen. Wie aus beiden Fig. 17 und 18 hervorgeht, 65 
ist es beispielsweise bei einer Transformationslange von N =» 512 giinstiger, bis zu 13 Spektralkoeffizienten uber 
den Goertzel-Algorithmus einzeln zu berechnen, als Qber eine schnelle Fourier-Transformation. Bei der Anwen- 
dung des Goertzel-Algorithmus in der BARK-Transformation gewinnt man daher einen ganz betrachtlichen 
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GeschwindigkeitsvorteiL Der Rechenaufwand fur das in Fig. 12 angegebene Beispiel reduziert sich dadurch 
enorm; in den meisten Fallen braucht sogar nur ein Spektralkoeffizient berechnet zu werden. Die Berechnung 
einer schnellen Fourier-Transformation uber eine Koeffizientenanzahl, die nicht einer Potenz von 2 entspricht, 
ist zwar m6glich, fuhrt aber im allgemeinen zu einem unverhaltnismaBig hdheren Aufwand aufgrund der 

s komplizierteren Struktur der dabei einzusetzenden Verfahren. 

Ein weiterer Vorteil des Goertzel- Algorithmus bestehtdemgegenaberdarin,daB beiiebige Frequenzrasterun- 
gen Uber dasselbe Verfahren (d h. uber das gleiche Programm) gerechnet werden kdnnen. Dazu ist lediglich die 
LSnge des Eingangsvektors und der Winkei 2nk/N gemaB der zu berechnenden Linie k neu vorzugeben. Obwohl 
der Goertzel-AIgorithmus a!s rekursives Filter am Stabilitatsrand arbeitet, treten dennoch keine Stabilitatspro- 

io bleme (z. B. Oberlauf des Rechenwerks) auf, da der Algorithmus nach der Berechnung eines jeden Koeffizienten 
der BARK-Transformation wieder neu aufgesetzt wird. 

FQr die BARK-Transformation ergibt sich somit gegenuber dem in Fig. 12 gegebenen Beispiel ein weitaus (urn 
den Faktor 224) geringerer Gesamtaufwand zur Ermittlung des gehorrichtigen Kurzzeitspektrums: 25 Goertzel- 
Algorithmen zur Berechnung der kompiexen BARK-Koeffizienten der 25 Teilbander und 1 1 Polyphasenfenster 

15 zur Einstellung der 1 1 verschiedenen Bandbreiten. 

Beide Teilalgorithmen erfordern im Kern ausschlieBlich Multiplikations-Additions-Operationen und sind 
daher far eine Ausfahrung auf einem Signalprozessor pradestiniert Die konsequente Anwendung des Prinzips 
der Reziprozitat und die Polyphasen-Filterung fuhrt fttr zu der in Fig. 20 fQr eine Abtastfrequenz von f a - 
44,1 kHz angegebenen Betragsubertragungsfunktionen der Einzeifilter einer BARK-Transformation, aufgetra- 

20 gen Ober einer linearen Frequenzachse normiert auf die Abtastfrequenz von f a - 44,1 kHz. Die DARK-Trans- 
formation weist ein konstantes Obertragungsverhalten im DurchlaBbereich, eine ausgezeichnete Sperrdamp- 
fung und eine hervorragende Weitabselektion (keine Nebenmaxima) auf. 

In Fig. 19 ist der Frequenzgang der Summe der Amplitudenubertragungsfunktionen angegebender BARK- 
Filterkaskade angegeben. 

25 

Besonderheiten der BARK-Transformation 

Die BARK-Transformation besitzt in der hier vorgestellten Form viele ideale Eigenschaften. Die Nachteile 
herk&mmlicher Frequenzzerlegungen uber Transformations- und Teilbandverfahren konnen damit uberwunden 
30 werden. Die Hauptvorteile der BARK-Transformation sind: 

- Gehorrichtige Frequenzauflosung: Die BARK-Transformation folgt dem von Zwicker ermittelten Zu- 
sammenhang far die Frequenzgruppenbreite an der jeweiligen Mittenfrequenz. Der Fehler, der durch die 
Vorgabe eines Frequenzrasters bei der iterativen Ermittlung der Grenzfrequenzen der frequenzgruppen- 

35 breiten Teilbander entsteht, kann beliebig klein gehalten werden, wenn nur das Raster fein genug ist 

- Feinstmdgliche Zeitaufldsung: Die Impulsantwort eines (Band-)Filters kann als Fensterfunktion far eine 
Transformation aufgefaBt werden. DieTransformationskoeffizientenstellen im Frequenzbereich Reprasen- 
tanten far die (Fourier-)Transformierten des Fensters im Zeitbereich dar. Deshalb muB sich bei der Varia- 
tion der "Bandbreite eines Koeffizienten" auch die Auflosung im Zeitbereich andern. Bei der BARK-Trans- 

40 formation ist daher far jede Frequenzgruppe das Produkt aus der Bandbreite B und der Zeitaufldsung T 
konstant und gleichzeitig minimal: B-T = 0,5. Diese Grenze, die als Heisenberg-Grenze bezeichnet wird, 
kann von keiner anderen Frequenzanalyse unterschritten werden. Der Filterentwurf entlang der Heisen- 
berg-Hyperbel gewahrleistet, daB weder die Zeitaufldsung noch die Bandbreite der Teilbander durch die 
BARK-Transformation verschlechtert wird. 

45 - Linearphasige Teilfilter: Aufgrund der Linearphasigkeit der Teilfilter und den gezeigten idealen Eigen- 
schaften der aus der Iteration resultierenden Betragsubertragungsfunktionen kann die BARK-Transforma- 
tion in dieser Form als bestens geeignet zur Nachbildung der Gehdreigenschaft angesehen werden, Fre- 
quenzgruppen zu bilden. Die Hnearphasigen Teilfilter erfordern keine Laufzeitentzerning durch Allpasse. 
Da aberdies die Filteriangen fur alle Teilfilter gieich sind, erubrigen sich auch Laufzeitanpassungen durch 

50 Verzftgerungsglieder. 

— Rechteckfdrmige Teilbander: Der iterative Filterentwurf aber die Fenstermethode kann interaktiv an 
einem Rechner erfolgen. Der Entwurf gewahrleistet - trotz der Beschneidung der Spaltfunktion auf eine 
endliche Lange — Betragsubertragungsfunktionen fur die frequenzgruppenbreiten Teilbander, die sich 
durch einen konstanten DurchlaBbereich und eine gute Weitabselektion wegen der fehlenden Nebenmaxi- 

55 ma auszeichnen. AuBerdem wird ein nahezu konstantes Obertragungsverhalten der gesamten Filterbank 
erreicht Zu den ausgezeichneten Eigenschaften, vor allem der hohen Flankensteilheit bei relativ geringem 
Realisierungsaufwand gegenaber einzelnen Bandfiltern, tr3gt nicht zuletzt die Polyphasen-Filterung bei. 

— Okonomische Realisierung: Bei der Reaiisierung wird ausgenutzt, daB die Frequenzgruppenbreiten im 
unteren Frequenzbereich nahezu konstant sind. Es muB daher mit weit weniger Fensterfunktionen im 

60 Zeitbereich multipliziert werden, als Teilbander nach der Transformation entstehen. 

Die mit der BARK-Transformation gefundene Frequenzanalyse aber nahezu rechteckfarmige Teilbander 
bildet den idealen Ausgangspunkt fur die weitere Entwicklung eines gehfirrichtigen MeBverfahrens einerseits 
und gehorangepaBter Codierverfahren andererseits. Das Transformationsverfahren laBt sich dariiber hinaus 
65 auch auf andere Anwendungen abertragen, bei denen einen Transformation aquidistanter Abtastwerte vom 
Zeitbereich in den Spektralbereich bendtigt wird, wobei die Spektralkoeffizienten auf der Frequenzachse 
nicht-flquidistant angeordnet sind und dennoch die Heisenberg-Bedingung erfOllt bleiben soil 
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Patentanspruche 

1. Frequenzanaiyseverfahren 

bei dem aus einem zeitdiskret in Form von Abtastwerten vorliegenden Eingangssignai mittels Fensterfunk- 
tionen uberlappend aufeinanderfolgende BIdcke von Abtastwerten entnommen werden, worauf die so mit 5 
den Fensterfunktionen bewerteten. BIdcke zur Erzeugung eines Satzes von Koeffizienten furier-transfor- 
miert werden, 
dadurch gekennzeichnet, 

daB jeweils ein solcher Block des Eingangssignals mit einem Satz von unterschiedlichen, symmetrischen 
Fensterfunktionen gleicher Lange bewertet wird, um einen Satz von mit den jeweiligen Fensterfunktionen 10 
gewichteten B16cken von Abtastwerten zu erzeugen, wobei sich die Fensterfunktionen durch die unter- 
schiedlichen Bandbreiten ihrer Fourier-Transformierten voneinander unterscheiden, und 
daB jeder der zeitgleich erzeugten BIdcke von gefensterten Abtastwerten eines Satzes einer eigenen 
Fourier-Transformation unterworfen wird, 

— indem zu jedem der zeitgleich erzeugten Blocke wenigstens ein Koeffizient berechnet wird, dem die 15 
Bandbreite der Fourier-Transformierten dieser Fensterfunktion zugeordnet ist, und 

— indem die Koeffizienten so gewahlt werden, daB die ihnen zugeordneten Frequenzbander sich, 
bezogen auf ihre Bandbreite, nur geringfiigig uberlappen oder nur geringfiigig voneinander beabstan- 
det sind. 

2. Frequenzanaiyseverfahren nach Anspruch 1, dadurch gekennzeichnet, daB die Fensterfunktionen sinx/ 20 
x-Funktionen sind, deren Mulldurchgange in Abtastzeitpunkte von Abtastwerten gelegt sind, so daB sich 
eine Zeitauflosung entsprechend dem gegenseitigen Abstand der Nulldurchgange ergibt, wobei die Zeitauf- 
losung umgekehrt proportional zur doppelten Bandbreite der Fourier-Transformierten ist 

3. Frequenzanaiyseverfahren nach Anspruch 2, dadurch gekennzeichnet, 

daB die Fourier-Transformation eine diskrete Schnelle Fourier-Transformation (FFT) ist, die eine der 25 
Anzahl von Eingangswerten im Zeitbereich entsprechende Anzahl von Koeffizienten im Frequenzbereich 
liefert, deren Abstand auf der Frequenzachse der Abtastfrequenz geteilt durch die Lange der Transforma- 
tion entspricht und deren Bandbreite durch die Fensterfunktion vorgegeben ist, und 

daB aus diesen Koeffizienten nur solche selektiert werden, die zusammen mit zeitparailel aus dem Satz der 
aus weiteren schnellen Fourier-Transformationen auf demselben Block von Abtastwerten mit unterschiedli- 30 
chen Fensterfunktionen zeitparailel zur Verfugung stehenden Koeffizienten bei anderen Mittenfrequenzen 
sich eine Kaskade von Bandfiltern aufbauen laBt, so daB die sich ergebenden Frequenzbander im wesentli- 
chen luckenlos aneinander angrenzen. 

4. Frequenzanaiyseverfahren nach Anspruch 2, dadurch gekennzeichnet, 

daB zur Verbesserung der Selektionseigenschaften der so entstehenden Kaskade von Bandfiltern die Lange 35 
der unterschiedlichen Fensterfunktionen eines Satzes auch ein Vielfaches der Transformationslange betra- 
gen kann, 

wobei vor der Transformation die Vielfachheit von Eingangswerten durch Aufteilung und additive Oberla- 
gerung (Superposition) der bereits gefensterten Eingangswerte wieder auf die Transformationslange redu- 
ziert wird (Polyphasenaddition). 40 

5. Frequenzanaiyseverfahren nach Anspruch 2, dadurch gekennzeichnet, daB die selektierten Koeffizienten 
aus dem Eingangssignai mittels eines Goertzel-Algorithmus zweiter Ordnung berechnet werden. 

6. Frequenzanaiyseverfahren nach einem der Anspruche 1 bis 5, gekennzeichnet durch 

ein iteratives Verfahren zur Bestimmung derdiskreten Mittenfrequenzen und Bandgrenzen zur Festlegung 
der Fensterfunktionen, welches folgende Schritte umfaBt: 45 

— Setzen einer Lauf frequenzvariabien (f m itte) auf eine Startfrequenz (foben). die auf dem durch die 
Lange der Fourier-Transformation gewahlten Frequenzraster liegt; 

— Verschieben der Lauf frequenzvariabien (fmitte) auf dem Frequenzraster um einen Schritt (df); 

— Ermittlung der sich aufgrund einer gegebenen, nicht linearen Funktion der Mittenfrequenz ergeben- 
den Bandbreite; und 50 

— Wiederholung der Iteration soiange, bis eine Bandgrenze des zu der Lauffrequenzvariablen (fmitte) 
gehorenden Bandes mit einer Bandgrenze eines vorher ermittelten Bandes zusammenfailt 

7. Frequenzanaiyseverfahren nach Anspruch 6, dadurch gekennzeichnet, daB die gegebene, nichtlineare 
Funktion der Zusammenhang zwischen der Mittenfrequenz und der Frequenzgruppenbreite ist 

55 
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Vorgaben: 
Abtastfrequenz: fa 
Frequenzraster: #lines 



T 




foben = fa / 2 
df = foben / #lines 



T 



fmitte = foben 



Tirir; 
T 



fmitte = fmitte - df 




funten = fmitte - (foben - fmitte) 



BW = (foben - funten) | 



foben = funten 




Ausgabe an Filterentwurf: 
BW(fmitte) 



I 
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Audio signal frequency analysis method • using window functions to provide sample signal 
blocks subjected to Fourier analysis to obtain respective coefficients. 
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The frequency analysis method involves dividing an input signal, represented by discrete sample 
values, into overlapping blocks using a window function and subjecting them to a Fourier 
transformation, to obtain a set of coefficients. Each input signal block is evaluated via a set of 
differing symmetrical window functions of equal length, their Fourier transformations having 
different bandwidths. 

The simultaneously generated blocks are used to provide a set of characteristic Fourier 
transformations, with a coefficient provided for each block, dependent on the window function 
Fourier transformation bandwidth. The selected coefficients have frequency bandswhich are only 
slightly overlapping, or only slightly spaced apart. 

USE/ AD VANTAGE - Accurate simulation of human hearing characteristics. At least first five 
spectral coefficients processed together, e.g. in fast fourier transformation. 

Dwg.6/20 

US 5583784 A 

A frequency analysis method that simulates perception characteristics of human hearing, comprising 
the steps of: 

extracting from an input signal a sequence of blocks of discrete sampled values, wherein successive 
blocks of the sequence overlap each other in time; 

concurrently windowing one of the successive blocks using a plurality of differing, 
symmetricalwindow functions of substantially identical length, each window function having a 
corresponding Fourier transform with a different bandwidth, to produce a plurality of windowed 
blocks of weighted sampled values; and 

performing a Fourier transformation on each windowed block wherein for each windowed block at 
least one coefficient is generated having assigned thereto the bandwidth of the Fourier transform of 
the window function used to produce the windowed block, and wherein the at least one coefficient for 
each windowed block is chosen such that bandwidths assigned to the chosen coefficientsof the 
plurality of windowed blocks adjoin one another approximately without gaps and in a manner 
approximately free from overlapping. 
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